Active diffusion of motor particles 
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The movement of motor particles consisting of one or several molecular motors bound to a cargo 
particle is studied theoretically. The particles move on patterns of immobilized filaments. Several 
patterns are described for which the motor particles undergo non-directed but enhanced diffusion. 
Depending on the walking distance of the particles and the mesh size of the patterns, the active 
diffusion coefficient exhibits three different regimes. For micrometer-sized motor particles in water, 
e.g., this diffusion coefficient can be enhanced by two orders of magnitude. 

PACS numbers: 87.16.Nn, 81.07.-b, 05.40.-a 



43 

o 



c3 



-a 
c 

o 
o 



> 
o 

(N 

in 
o 

-I— > 
c3 



i 

G 

O 

o 



X 



Introduction. Biological motor molecules such as the 
cytoskeletal motors kinesin and myosin, which convert 
the chemical free energy released from the hydrolysis of 
adenosine triphosphate (ATP) into directed movements 
along cytoskeletal filaments, represent powerful nanoma- 
chines. It seems rather promising to construct artificial 
biomimetic devices based on these molecular motors and 
to use them, e.g., for the transport of cargo particles 
in the nano- and micrometer regime Q, Q- These ap- 
proaches are based on two different in vitro motility as- 
says, the gliding assay, for which molecular motors are 
anchored at a substrate surface and drag filaments along 
this surface, and the bead assay, for which the filaments 
are immobilized on the surface and the motors pull beads 
or cargo along these filaments Q. 

In order to obtain spatial control of the movements, 
motor molecules or filaments can be bound to the sub- 
strate surface in well defined patterns by controlling 
the surface topography and/or chemistry. Examples are 
motor-covered surface channels or grooves which guide 
the gliding filaments, chemical patterning of a surface 
to bind filaments or motors selectively, and combina- 
tions of topographic and chemical patterning, see, e.g., 
[H @ S Etl In this way, one can create 'ac- 
tive stripes' along which the active transport takes place. 
The directionality of these stripes can be controlled in a 
rather direct way for systems based on the bead assay, 
since unidirectional movement of motors can be realized 
by alig ning filaments in parallel with the same orientation 
yjllfj. In gliding assays, movement of filaments usually 
occurs randomly in both directions, but unidirectional 
movements can be obtained constructing special surface 
domain patterns Q. 

Usually, molecular motors undergo directed movements 
which can be used for transport over large distances. 
However, the active motion of molecular motors can also 
be used to generate effectively diffusive movements by 
sequences of active directed movements into random di- 
rections. We will call these ATP-dependent and thus 
energy-consuming, but effectively non-directed move- 
ments active diffusion. It should be useful under con- 
ditions where simple diffusion is slow, i.e., for large car- 



goes or high viscosity of the fluid. Active diffusion is also 
present in vivo, where various motor-driven cargoes such 
as vesicles, messenger RNA, and viruses exhibit active 
back-and- forth movements . 

In the following, we will discuss active diffusion in sev- 
eral simple systems which can be realized experimentally. 
We present results for several cases where patterns of 
'active stripes' have been created on a surface by im- 
mobilizing cytoskeletal filaments in certain well-defined 
patterns, so that they constitute a system of tracks for 
the active movements of particles coated with molecular 
motors. Even though we focus on these latter systems 
derived from bead assays with mobile motors, our argu- 
ments can be easily extended to mobile filaments gliding 
in motor-covered surface grooves. 

Model. In order to study the effective diffusion aris- 
ing from active motion along patterns of filaments, we 
describe these movements as random walks on a lattice 
[lil lull ]. We consider a cubic lattice with lattice con- 
stant I on which certain lines of lattice sites represent 
the filaments or active stripes. If a motor-driven parti- 
cle occupies such a lattice site, it makes a forward step 
(along the orientation of the filament), a backward step, 
and no step with probabilities a, (3, and 7 per time unit 
t, respectively. In addition, it may escape to each of 
the n a( j adjacent non-filament sites with probability e/6 
("■ad = 3 and 4 for a filament at a wall and in solution, 
respectively). An unbound particle, i.e., a particle at 
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FIG. 1: Lattice model for the active motor movements along 
filaments. Motors make forward and backward steps along 
the filament with probabilities a and /3, respectively and no 
step with probability 7. Unbinding from the filament occurs 
with probability eo = Wadt/6 where n a( j is the number of non- 
filament neighbor sites, and a motor reaching the filament 
binds to it with sticking probability 7r ad . 
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FIG. 2: Active diffusion in a tube with two antiparallel fil- 
aments (double- logarithmic plot): The variance Ax 2 of the 
particle position as a fuction of time t for motor-driven par- 
ticles with unbinding probabilities e = 0.01 (solid), 10~ 3 
(dashed), and 10 -4 (dotted line) and other parameters as de- 
scribed in the text. The thin dashed line indicates unbound 
diffusion which is much smaller than active diffusion. Here 
and in the following figures, lengths and times are expressed 
in units of I and r, respectively. 



a non-filament site, performs a symmetric random walk 
and steps to each neighbor site with probability 1/6 per 
time t. An unbound motor reaching a filament site binds 
to it with probability 7r ac i . The parameters of the random 
walks can be adapted to the measured transport proper- 
ties 0] ■ Confining surfaces are implemented by rejecting 
all stepping attempts to lattice sites of the surfaces. 

For the simulations, we chose the basic length scale to 
be I = 50nm and the basic time scale to be r = 25ms. 
These values result in a diffusion coefficient £> u b ~ 
10 _1 /im/s as appropriate for a micrometer-sized parti- 
cle in water. We take a — 0.5, (3 — 0, and 7r ac i = 1, 
corresponding to a velocity t>b — 1/im/s and a binding 
rate in the order of 100s _1 for motors within a capture 
range of 50nm from the filament. 

Active diffusion in Id. In order to discuss the simplest 
case of active diffusion, let us consider two parallel fila- 
ments with opposite polarity. We denote the coordinate 
axis parallel to the filaments by x and take the filaments 
to be sufficiently long, so that the filament ends are not 
reached within the observation time. The filaments are 
placed within a long and thin tube or channel with a 
rectangular cross section (with width L and height H), 
so that unbound motors cannot escape, but rebind with 
a finite mean return time, a situation that mimics the 
axon of a nerve cell. 

We have studied this case u sing an analytical Fourier- 
Laplace technique as in Refs. [13| and computer simula- 
tions. Some results are shown in Fig. [3 Because of the 
antiparallel orientation of the filaments, no net movement 



is obtained in such a system, and the effective velocity 
is zero on large scales. The effective diffusion coefficient, 
D = lim^oo Ax 2 (t)/2t, defined via the long-time be- 
haviour of the positional variance Ax 2 = (x 2 ) — (x) 2 , is 
found to be given by 

D = D + D act P h , (1) 

where the first term is the weighted average of the bound 
and unbound diffusion coefficient, D — D^P^ + Z? u b(l — 
P h ) with the probability P h = 1/[1 + {e/n^LH/e 2 ] that 
the motor is bound to the filament at large times. The 
contribution of the bound diffusion coefficient is usually 
small, so that D ~ Z? u b(l — Pb)- The second term de- 
scribes active diffusion, i.e., additional spreading of a 
distribution of many motors due to the active directed 
movements into both directions with 

1v 2 

= ^. (2) 

The latter contribution can be written as D act = 
i(Lg)/Ai b = Ax\yV h with the average walking time 
At h = 3/(2e), the average walking distance A^b = 
fbA^b, and the mean square of the distance L s traveled 
actively before a change of direction via unbinding, as 
obtained from the distribution P(L S ) — 2K^~ e ~^ Ls ^ AXb 
of the walking distances. 

This contribution is present even if the diffusive motion 
of unbound motors is strongly suppressed, e.g. because 
of their large cargoes or because of a high viscosity of 
the fluid. In cells, this should apply to the movements of 
large organelles such as mitochondria (where, however, 
the switching of direction is accomplished via different 
types of motors rather than via the orientation of the 
filaments) 0]. 

Active diffusion as given by Ax 2 w 2Dt is valid at large 
times with t > t* = 2D/i% ~ 3P h /e, see Fig. H Note 
that the larger the active diffusion coefficient, the larger 
is also the crossover time i* and, thus, the time necessary 
to observe the active diffusion. 

An enhancement of diffusion occurs also in rather dif- 
ferent contexts such as electrophoresis, chromatographic 
columns, and diffusion in hydrodynamic flow which can 
also be described by multi-state random walk models 
0,0]. In all of these cases, however, the directed mo- 
tion arises from externally applied fields rather than from 
the active movements of the motor particles. 

Arrays of filaments on a surface. Let us now consider 
filaments which are immobilized on patterned surfaces 
with two-dimensional arrays of stripes. We consider four 
systems as shown in Fig. |3J (Al, A2) consist of arrays 
of parallel stripes whereas (Bl, B2) correspond to in- 
tersecting stripes which form a two-dimensional square 
lattice with lattice constant L. The stripes have a width 
W; in our simulations, we used W = I. In both cases, 
we studied two situations: (i) Each stripe contains fil- 
aments with both possible orientations as in (Al, Bl), 
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and motor particles which bind to a stripe can start to 
move in either direction (which is chosen randomly in the 
simulations); and (ii) neighboring stripes are covered by 
filaments with opposite orientation as in (A2, B2). Note 
that all patterns have no directional bias on large scales. 

The motors can change their direction of movement 
via unbinding and rebinding to another filament. In the 
stripe lattices, they can also change direction at the ver- 
tices by switching to another filament, which is chosen 
randomly. Finally, we take the diffusion in the direction 
perpendicular to the surface to be restricted to a slab of 
height H . To study the active diffusion for this case, we 
determined the diffusion coefficient by extensive Monte 
Carlo simulations for motor-driven active particles of mi- 
crometer size and for a wide range of the unbinding rate 
e and the stripe separation L as shown in Fig. 0] 

Arrangements of filaments of this type are accessible to 
experiments and can be created using structured surfaces 
as described above. Lattice patterns of filaments can also 
be formed by the self-organization of filaments and motor 
complexes [l6j or by the assembly of crosslinked networks 
on micropillars 17]. 

Parallel active stripes. Arrays of parallel active 
stripes which contain filaments with both orientations 
(Al) exhibit active diffusion very similar to the one- 
dimensional case discussed above. However, in this case, 
the effective diffusion coefficient depends on the spatial 
direction. While the effective diffusion is described by 
Eq. in the direction parallel to the stripes, which leads 
to the parallel diffusion coefficient Dm = D + D act Pb, the 
diffusion perpendicular to the stripes is characterized by 
D± = D u b(l — Pb), see Fig. ^a). The latter expression 
is also valid if the directionality of the stripes alternates 
(A2), where the unbound motors return predominantly 
to the same stripe which leads to a higher value of Dm as 
shown in Fig.^a). 

Regular lattices of active stripes. In the system 
architectures (B1,B2), the motor-driven active parti- 
cles can perform two-dimensional random walks on the 
stripe lattice even if they remain bound to the fila- 
ments. In addition, effectively diffusive behavior can 
arise from the repeated unbinding from and binding to 
filament. Let us consider case (Bf). The diffusion co- 
efficient is given by D — D act Pb + D u b(f — Pb) with 
P b « l/[l + {e/ix &A )LH/{2Wl)] for the stripe lattice. De- 
pending on the processivity of the active particles, i.e. on 
the value of the unbinding rate, there are three regimes: 

(I) If the unbinding rate is small, so that the walk- 
ing distance is large compared to the pattern mesh size, 
Axb = 2i>b/e 3> L, the particle performs a random walk 
along the active stripe lattice and the diffusion coefficient 
is D w D(e = 0) = jv h L. This case should be appropri- 
ate for cargo particles driven by many processive motors. 

(II) If the average walking distance before unbinding 
is smaller than the mesh size, Axb = 2t>b/e < L as ap- 
propriate, e.g., for transport by single processive motors, 
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FIG. 3: Substrate surfaces with active stripes which consist of 
filaments immobilized onto striped surface domains: For the 
system architectures (AI) and (Bf), all stripes contain fila- 
ments with both orientations, for (A2) and (B2), the filaments 
within one stripe have the same orientation, but filaments on 
neighboring stripes have opposite orientation. 



changes in direction of movement will mainly be through 
unbinding and rebinding to another filament. In that 
case, the distance traveled before a change of direction 
is given by the walking distance A^b and, similar to the 
one-dimensional case discussed before, the diffusion coef- 
ficient is D w iAxbWbPb 0- 

(III) Finally, if motor particles spend most of the time 
unbound, i.e., for large e as appropriate for weakly pro- 
cessive motors, the diffusion coefficient is simply given by 
the unbound Brownian motion, D « D u b- 

The three regimes can be summarized in one formula 
by using D act = VbL s /4 and defining the interpolated 
distance traveled between two changes of direction by 
L s = L2Axb/(L + 2Axb) which leads to 



D = 



LAxv, 



2 L + 2Ax] 



-Pb + D ub (f-P b ). 



(3) 



The latter expression leads to the solid line in Fig.^Jb), 
which agrees well with the data points (open circles) from 
our simulations. 

The value of the active diffusion coefficient can be fur- 
ther enhanced by controlling the orientation of the fila- 
ments along the stripes, i.e., by using filament patterns 
such as (B2) for which we obtain an active diffusion co- 
efficient twice as large as for case (Bf). In regime (I) 
this factor is due to the fact that a fully random choice 
of the new direction is completed after a walk over the 
distance 2L. However, Eq. with an additional factor 
of two in the active diffusion term describes rather well 
our simulation data also in regime (If). 

Discussion. Finally, let us estimate the order of 
magnitude of the active diffusion in systems that are 
experimentally accessible. Free diffusion in a Newto- 
nian fluid is governed by the Stokes-Einstein relation 
from which the diffusion coefficient is obtained as D u b — 
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FIG. 4: (a) Effective diffusion coefficients Du (circles) and D± 
(triangles) as functions of the stripe separation L. The open 
and filled symbols correspond to the two system architectures 
(Al) and (A2), respectively, compare Fig. [3] The full lines 
indicate the analytical result for (Al). (b) Effective diffusion 
coefficient D as a function of the unbinding probability e. 
D exhibits three regimes (I)— (III) as discussed in the text. 
The open and filled symbols correspond to the two system 
architectures (Bl) and (B2), respectively, as in Fig. [3] The 
parameters are as described in the text; we used e = 0.0025, 
H = 100£ in (a) and L = H = lOOi in (b). 

(100nm/i?hyd) x 2.4/im 2 /s x ry W atcr/?7- In the latter ex- 
pression, i?hyd is the effective hydrodynamic radius of 
the particle and 77 the viscosity of the medium. For 
a micrometer-sized particle in water, we have Z? u b ~ 
0.1/im 2 /s, for a lOOnm-sized particle D u b ~ l/im 2 /s. 
(The diffusion coefficient of bound motors moving along 
filaments is much smaller, typically of the order of Z?b ~ 
Vht ~ 0.01/im 2 /s with Vb ~ 1/zm/s and I ~ lOnm.) 

The active diffusion coefficient, on the other hand, is 
of the order of jV\>L s Pb with L s given by the filament 
length in regime (I) and by the walking distance in regime 
(II). With a typical filament length of a few tens of /zm, 
this means that active diffusion coefficients of the order 
of 10/im 2 /s can be obtained for cargo particles driven by 
many motors, for which regime (I) is appropriate. For 
micrometer-sized particles this is two orders of magnitude 
larger than the diffusion coefficient of unbound Brownian 
motion. For cargo particles driven by a single or a few 
motors, regime (II) applies, and for a typical unbinding 
rate of 1/s, we obtain an active diffusion coefficient of 
~ l/zm 2 /s. 

These estimates as well as our simulation data are for 
diffusion in water. In a more viscous medium with viscos- 
ity 77 > f7 W atcr , diffusion is even more enhanced due to the 
active processes, because active diffusion does not obey 
the Einstein relation and is largely independent of 77 |l9| : 
the viscous force on a micrometer-sized bead moving with 
l/jm/s is of the order of ~ rj/rj watcv x 10 _2 pN which is 
small compared to the piconewton forces required to slow 
down the active movements of a single motor. Unbound 
diffusion is however reduced by the factor 77 wat0 r/?7, so 
that the ratio D/D u b ~ 77 increases with increasing vis- 
cosity 77. 

Active diffusion will therefore be useful for particles 



which explore surface patterns with linear dimensions of 
tens of micrometers or even millimeters, e.g., in order to 
find a specific binding partner or to deliver a cargo, in 
particular if the surrounding fluid has a high viscosity. 
One potential application is provided by the integration 
of these biomimetic transport systems into labs-on-a-chip 
for genomics or proteomics. 

In summary, we propose here to use chemically or to- 
pographically patterned surfaces to create well defined 
arrangements of active stripes covered by filaments and 
to use the random walks arising from the active move- 
ments of motor-covered particles along these stripes to 
achieve fast diffusion on the surface. 
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